Finite temperature spectral function of Mott insulators and CDW States 
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We calculate the low temperature spectral function of one-dimensional incommensurate charge 
density wave (CDW) states and half-filled Mott insulators (MI). At T = there are two dispersing 
features associated with the spin and charge degrees of freedom respectively. We show that already 
at very low temperatures (compared to the gap) one of these features gets severely damped. We 
comment on implications of this result for photoemission experiments. 

PACS numbers: 71.10.Pm, 72.80.Sk 



The determination of the finite temperature single- 
electron Green's function in strongly interacting one- 
dimensional (ID) electron systems with a spectral gap 
is a problem of significant interest as the spectral func- 
tion A(u,q) is expected to exhibit the celebrated prop- 
erty of spin-charge separation. Rather than a single co- 
herent quasiparticle peak, one expects to see two broad 
features in A(u>,q), which are associated with the inde- 
pendent dispersion of spin and charge collective modes. 
This phenomenon also occurs in the case of a ID metal- 
lic Luttingcr liquid state. The great advantage of the 
insulating state is that due to the presence of a gap it 
is much more robust against the effects of 3D couplings, 
temperature or impurities. Hence the chances to observe 
spin-charge separation by means of Angle Resolved Pho- 
toemission Spectroscopy (ARPES) are higher in Mott or 
CDW insulators than in metallic systems. Unfortunately 
it is much more difficult to determine dynamical correla- 
tion functions in these states. The very existence of a gap 
prohibits the application of methods based on conformal 
field theory such as the Luttingcr liquid approach. A 
further complication is that the gap is dynamically gen- 
erated, which invalidates mean-field approximations to 
the problem. In light of these difficulties, 1/N-expansions 
have been a method of choice, where A is a large parame- 
ter introduced by enlarging the spin rotational symmetry 
of the Hamiltonian from SU(2) to SU(N). 
In this letter we derive asymptotically exact expressions 
for the spectral function and the tunneling density of 
states for incommensurate CDW states and half-filled MI 
in the field theory limit. The reason to treat both cases 
at once is that the results for the MI can be obtained 
from those for the CDW state. At T = we perform 
our calculations for general N, which allows us to see 
how the large- A result progressively loses its accuracy 
as N decreases and for N = 2 no longer reproduces even 
qualitative features of the solution. We demonstrate that 
Luttinger's theorem continues to hold despite the absence 
of a Fermi surface. Finally we consider the effects of 
temperature in the physically relevant case N = 2 by 
means of a systematic expansion in exp(— A/T), where 
A is the spectral gap. There has been much previous 
work on determining the T = spectral function of ID 



MI and CDW states. In Ref a conjecture for A(u>,q) 
was put forward, which agrees with the result we obtain 
in the field theory limit by means of an exact, system- 
atic method. Expressions for the spectral function have 
also been obtained in the limit where the single-particle 
gap is much larger than the bandwodth 0. This regime 
is complementary to the case we address here. Finally 
there also have been extensive numerical studies on t-J 
and Hubbard models e.g. Refs [J 9- 
Incommensurate CDW state: As a starting point for a 
microscopic description of the CDW state one may choose 
a model of noninteracting electrons at some incommen- 
surate band filling, coupled to ID phonons. Examples 
are the Su-Schrieffer-Heeger and Holstein Q Hamil- 
tonians. The electronic low-energy degrees of freedom 
have momenta close to the Fermi momenta ±kp. In the 
low-energy sector the electron operators can therefore be 
represented as 



ikpx 



R a (x) + e 



-ikpx 



L a [x) 



(i) 



where do is the lattice spacing, x = nao and R and L are 
slowly varying Fermi fields. In three spatial dimensions 
the presence of an electronic spectral gap would auto- 
matically imply the formation of an anomalous average 
(R\,(x)L (J (x)) ^ 0. In one dimension the average is not 
formed even at T — 0; instead correlation functions of 
the operator R\L a decay in a power-law fashion. In the 
limit of infinite phonon frequency luq the phonons can 
be integrated out without inducing retardation effects in 
the resulting effective electron Hamiltonian 

7]- 

The case 

of large but finite ojq can be treated similarly, as long 
as one is only interested in low energies w C uo i- In 
this regime the electronic model obtained by integrating 
over the phonon degrees of freedom is described by the 
following universal Hamiltonian 



H = Hc + Hs 



(2) 



H c = ^[K-^d^f + K^dM 2 ], (3) 

I D7T 



H s = 



[: J a J a : + : J a J a :} + v s gJ a J a . (4) 



Here $ c is a canonical Bose field, C is its dual field and 
J a = L\r% a iL a >, J a = R) a T^ a ,R a i are current operators 



2 



satisfying the level- 1 SU(2) Kac-Moody algebra. The pa- 
rameters v c .s (charge and spin velocities), A (spin gap) 
and K c (Luttinger liquid parameter) depend on the de- 
tails of the underlying microscopic lattice model and in 
particular on uig. The parameter K c controls the scaling 
dimensions in the charge sector. 

Half-filled MI: The Hamiltonian thirring-curr is identi- 
cal to the low-energy theory for a half-filled MI like the 
Hubbard model, provided we interchange charge and spin 
sectors c <-*■ s and then set K s = 1, kp = tt/2. Below we 
present results for the more general CDW case with the 
understanding that the corresponding results for the MI 
are obtained by the aforementioned mapping. 
Zero Temperature: The model thirring-curr is general- 
ized to SU(N) by replacing the currents in curr by their 
SU(N) analogs. The spectrum in the charge sector is 
gapless E c (k) — v c \k\, whereas the spectrum in the spin 
sector consists of scattering states of electrically neutral, 
massive solitons with dispersion E s (k) = -\/A 2 + (v s k) 2 
and their bound states for N > 2; the presence of bound 
states does not affect the results presented here. For 
small g <C 1 we have A = Dg 1 / N e~ 2n / Na , where D ~ ujq 
is the ultraviolet cut-off. The electron operators can be 
expressed as products of (vertex) operators acting in the 
charge and spin sectors respectively 

(L a (r,x) 4(0)) = U(0 (r l8 ) O*(0)). (5) 

Q — C,S 

Using the results of [lfllll] ] we obtain the following result 
for the asymptotics of the single-particle Green's function 
at zero temperature 



G(r,x) = 


e- ikp "G H (r,x)+e ikFX 


G L ( 


T,x) , 


(6) 


G l (t,x) ~ 


Z N A 


2v c /A 


1 
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[(2«c/A) 2 " 
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y/TTV s VcZTr 


{v c t + ix) 




X 2 + V 2 T 2 




X 


( v s t — ix\ 


V — 1 
2N 


(Ar) , 


(7) 



V S T + IX 



where r = \/ t 2 + x 2 v s 2 
K, 



and & = jn 



J- - 2 



Zjy is a nonuniversal constant 
We note that Gr(t,x) = 

Gl{t, —x). Equation TO is derived by taking into ac- 
count processes with emission of an arbitrary number 
of charge excitations and only one massive spin excita- 
tion. The corrections to TO are of order 0(e~ 3Ar ) for 
N = 2. However, due to the fact that matrix elements 
corresponding to the emission of multi-soliton states are 
numerically small, the accuracy achieved by this approx- 
imation is very good. The smallness of matrix elements 
also allows us to generalize our calculations to finite T. 
In the case of equal velocities v s = v c = v the retarded 
Green's function at T = for right moving fermions is 



uj + vq f 9 1 8 s > 
X V^( 1 -2' 2 -iV-2< 2; A^' (8) 

where F is a hypergeometric function and s 2 = lu 2 — 
(vq) 2 . As might be expected, the gap in the spectral func- 
tion is 'clean', i.e. Ar(oj, kp + q) = -ImGj{(w, kp + q)/ir 
vanishes for \ui\ < y A 2 + v 2 q 2 . The zero temperature 
tunneling density of states is 



A 



N 



acosh(o;/ A) 



dx 



cosh(x[l -1/N}) 
(w/A-coshx) 1 - 9 - 1 /^ 



(9) 



where A N = ^v s vJxt>+i/N ) We see that at T = P( u ) 
vanishes inside the gap. The singularity just above the 
threshold (0 < (w/A) - 1< 1) is 

A N B[ 



N ' 2/ 



V2 



(w/A- l) e+ ^ 



(10) 



We note that for the physical case N = 2 the "asymp- 
totic" result power is actually equal to fhoN. According 
to 11 0|l the behavior of p{ui) above the gap is determined 
by the scaling exponent 1/2 — 1/N+8. The non-universal 
part 8 is small for large phonon frequencies. On the 
other hand, the remaining part is determined only by 
the Lorentz spin of the spinon creation operator, and 
therefore is universal. Equations fhoN-power show that 
a 1/iV-expansion fails completely in the physically rele- 
vant case N = 2, where the tunneling density of states 
experiences a qualitative change and becomes nonsingu- 
lar at the threshold. As we have alluded to earlier, the 
TV — > oo limit agrees with the mean-field results of Ref. 

Luttinger's Theorem: The Green's function GrN has 
branch cuts, but no poles. In particular there are no 
poles at zero frequency and thus no Fermi surface. Nev- 
ertheless Luttinger's theorem as stated in Ref. |l2( is ful- 
filled since the logarithm of the Green's function In G(lu — 
0, kp+p) is still singular at the non- interacting Fermi sur- 
face, because the Green's function GrN has zeroes. This 
property follows from (i) the fact that the charge sector 
is gapless, which implies that (R a (T,x)Ll.(0)) = and 
(ii) Lorentz invariance of the low-energy effective theory, 
which implies that 

(* a (r,a;)*t(0)> ~ exp(±#)ft(r); $ = R,L. (11) 

Here r and <j> are polar coordinates and 1Z denotes the 
radial part of the correlation function. As we are dealing 
with an insulating state we have TZ(r) oc exp(— Ar) at 
large distances and hence J dr lZ(r)r is finite. Thus 



G«,i(0,0) = 



exp(±i</») / dr K(r)r = 0. (12) 



G R (w, kp + q) = —^T 



i-_]ri2-i-^ 

2 \ N 2 



For a metallic state the r integral would diverge and the 
Green's function would have a singularity rather than a 
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zero. 



Finite Temperature: Let us now turn to the general case 
v c 7^ v s and finite temperatures. We will restrict the dis- 
cussion to N — 2. The spectral function is conveniently 
expressed as a convolution 



Ar(u, k F + q) 



1 



(27T) 



dui'dq' g s (J ' ,q l ) 



x [g c {w -Lu',q- q') + g c {-w - lu' , -q - q')} , (13) 

where g CjS (w, q) are the Fourier transforms of the finite 
temperature correlators (0 CtS (x, t) O| s (0, 0))t- In the 
charge sector we can use the standard conformal mapping 
to obtain (see e.g. fl3j |) 



/(a, 7) 



C(T) f 



lu + v c q 1 + 1 
2nT 



Re 



{-2i)~<B 



2 

ia 



f 



v c q 



2vrT ' 2 



: ,l-7 



(14) 



where C{T) = (w c /27r 3 ) 1 / 2 (2 7 r/A) e T 6 ^ 2 and B(x,y) 
is the Euler Beta function. The correlation function in 
the spin sector can be determined by using exact results 
for the SU(2) Thirring model. We invoke a spectral re- 
spresentation in terms of scattering states of solitons and 
antisolitons, constructed by means of the Zamolodchikov- 
Faddeev algebra (see e.g. p. 7 of Refill. We intro- 
duce an index e = ± for solitons and antisolitons re- 
spectively and parametrize energy and momentum by a 
rapidity variable 9 in the usual way as E(9) = A cosh 9, 
P(0) = (A/v s ) sinh 6. In a basis of scattering states 
\6 n . . . ^i)e„... £l of solitons and antisolitons with rapidi- 
ties Oj the following spectral representation for thermal 
two-point functions holds 



n A< L j} i— 1 m d e ' fe } 

T~L TTL Tt 771/ 



x e 



fc=i 



fc=i 



As was shown recently, the representation specrep is 
suitable for carrying out a low-temperature expansion 
(T < A) of correlation functions fl5t fl6j|. Taking into 
account the two leading terms n = 0, m = 1 and n = 1, 
m — in specrep and combining this result with gc and 
all we obtain 



A R (uj,k F + q ) w A dz e z/2 



g c {uj - c(z),q-s(z)) 



+ e- c ^^ T g c (u + c(z) iq + S (z)) + 







f] 


u- 


* -9 J 





(16) 



3, 

EE 

< 




where A 



A cosh z and s(z) 



(A/v s ) sinhz. Let us now see how changing 9 and T af- 
fects An{io,q). We constrain our discussion to the cases 
of zero temperature and varying 9 and of finite T and 
fixed 9. As a function of 9, Au(ui,kF + q) displays a 
strongly varying behaviour. For small 9, correspond- 
ing to a high phonon frequency, the spectral function 
is very similar to the one of the half-filled MI (see the 
dotted curves in Fig. 2): there are two sharp, dispersing 
features associated with the spin and charge degrees of 
freedom respectively. For smaller phonon frequencies and 
concomitantly larger 9 these features become less promi- 
nent until they eventually disappear altogether. We plot 



CD/ A 



FIG. 1: Spectral function A R (u, k F + q) at T = 0, 6 = 0.8, 
v c = QAv s for v s q/ A = —2, —1,0, 1,2 (from bottom to top). 
The curves for different values of q have been offset. 



Ar(u), k F +q) at T = for the intermediate value 9 = 0.8, 
v c = 0.4 v s and several values of q in FigQ] We see 
that the peak associated with v s is already quite weak. 
The effects of a small, finite temperature T — 0.05A are 
shown in Fig[3for a half-filled MI (we now switch spin 
and charge sectors as discussed above and set 9 = 0, 
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kp = 7r/2). Compared to the T = result we find that 
the holon peak is significantly damped although T is still 
quite small compared to the low energy scale A. The 
physical reason for this is very simple: in the MI only the 
charge sector is protected by the gap, whereas the gapless 
spin sector is significantly affected by T. This leads to a 
damping of the peak associated with the charge degrees 
of freedom, whereas the spinon peak stays rather sharp. 



+ 

3 




where D = (2Tr/A) e (8n/v c ) 1 / 2 T 9 - 1 / 2 . In the regime 
T « w < A we may use contour techniques to extract 
the leading contribution to fhoT2 



Z,2" 



T 



-(A-«)/T 



(19) 



This shows that at low temperatures only an exponen- 
tially small fraction of spectral weight gets transferred 
into the gap. In summary, we have calculated the low 
temperature spectral function for incommensurate CDW 
states and half-filled Mott insulators in the field theory 
limit. Luttinger's theorem is shown to hold despite the 
absence of a Fermi surface. We have studied the ef- 
fects of temperature on the two dispersing features as- 
sociated with spin and charge degrees of freedom. We 
demonstrated that already a small temperature essen- 
tially wipes out the holon peak in the half-filled Mott 
insulator. We thank A. Chubukov, A. Finkelstein, P.D. 
Johnson and T. Valla for important discussions. This 
work was supported by the DOE under contract number 
DE-AC02-98 CH 10886. 
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FIG. 2: Spectral function for a half-filled MI (0 = 0, 
v s = 0.5v c ), T = (dotted) and T = 0.05A (solid) v c q/A = 
— 1.5, —1, . . . 2.5 (from bottom to top). The curves for differ- 
ent values of q have been offset. 

We note that ARPES measurements on ID cuprate 
Mott insulators are taken at room temperature in or- 
der to avoid charging of the sample. The temperature 
used in Figgis chosen to reflect the ratio T/A in the ex- 
periments on Sr2CuC-3. Although the theory presented 
here does not describe the ID cuprate MI quantitatively 
(U/t is too large in these compounds), we believe that 
our result gives a strong indication that temperature ef- 
fects are important and at least partially account for the 
fact that the holon feature is barely visible in the data. 
The low-T behavior of the tunneling density of states can 
be analyzed analogously. We find 



p{uj) 



dx 

-co (2^ 



p s (x)p c (oj - x) +lo —> —uj, (17) 



where 

p c {x) = -DRe 

p s (x) ~ Z 2 



{-2i) e+1 iB 



1 + 29 



uj 1 



2ttT' 2 



M«>-A) 



-M/T 



A) 



H-Vi 



vv; 



\uj\ + Vw 2 - A 2 

~Tu 



A 2 



. (18) 
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